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Abstract 

The group classification of models of axion electrodynamics with arbitrary 
self interaction of axionic field is carried out. It is shown that the extension 
of the basic Poincare invariance is possible only for constant and exponential 
self interaction terms. With using the three-dimensional subalgebras of the 
. Lie algebra of Poincare group an extended class of exact solutions for the elec- 

tromagnetic and axionic fields is obtained. These solutions include arbitrary 
parameters and arbitrary functions as well, the most general solutions include 
six arbitrary functions. Among them there are bound and square integrable 
solutions which propagate faster than light. 
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1 Introduction 



Group analysis of partial differential equations is a well developed branch of math- 
ematics including a number of interesting fundamental problems. But maybe the 
main its value are the powerful tools for construction of exact solutions of compli- 
cated nonlinear equations. Sometimes it is the group analysis which gives the only 
hope to obtain at least some solutions for complicated physical (chemical, biological...) 
models. 

In this preprint we present some results obtained with application of the group 
analysis to the important physical model called axion electrodynamics. There was a 
lot of motivations for this research, both physical and mathematical. 

Axions are hypotetical particles which belong to main candidates to form the dark 
matter, see, e.g. []]. Additional arguments to investigate axionic theories appeared in 
solid states physics, since the axionic-type interaction terms appears in the theoreti- 
cal description of crystalline solids called topological insulators [2] . The existence of 
a pseudoscalar (axion) field can be extracted from the experimental data concerning 
electric field induced magnetization on C^Os crystals or the magnetic field-induced 
polarization [3]. In addition, the axion hypothesis makes it possible to resolve the 
fundamental problem of quantum chromodynamics connected with the CP symme- 
try violation in interquark interactions [I]-[S]- Thus it is interesting to make group 
analysis of axionic theories which are requested in the three fundamental branches of 
physics, i.e., the cosmology, condensed matter physics and quantum chromodynamics. 

Let us present two more motivation with are very inspiring for us. Recently a new 
exactly solvable model for neutral Dirac fermions had been found and other integrable 
models for such particles had been indicated [5] . But these models involve an external 
EM field which do not solve Maxwell equations with physically reasonable currents. 
However, these fields solve equations of axion electrodynamics. 

We had classified exactly solvable quantum mechanical systems including shape 
invariant matrix potentials [9J, [ID] . Some of these systems also include solutions of 
field equations of axion electrodynamics. 

We have described the finite dimensional indecomposable vector representations 
of the homogeneous Galilei group and construct Lagrangians which admit the corre- 
sponding symmetries [H]-[l3] . And exactly the Lagrangian of axion electrodynamics 
appears to be the relativistic counterpart of some these Galilei invariant Lagrangians 

M- 

In addition, axion electrodynamics is a nice and complicated mathematical model 
which needs good group-theoretical grounds. In this preprint we present such grounds 
and also find an extended class of exact solutions for the related field equations. 



1 



2 Equations of axion electrodynamics 



Let us consider the following generalized Lagrangian: 

L = - + \QF^ - V(9). (1) 

Here F^ u = d^A u — d u A^, is the vector-potential of electromagnetic field, F^ u = 
\t[j,u P tjF pa ', 9 is the axion field, = <9 M #, V{9) is a function of 9 and k is a dimen- 
sionless constant which is supposed to be nonzero and can be rescaled to the unity. 

If 9 = then formula ([I]) gives the Lagrangian of Maxwell field. For V{9) = ^m 2 9 2 
equation ([T]) reduces to the standard Lagrangian of axion electrodynamics. 

Our goal is to make the group classification of the Euler-Lagrange equations cor- 
responding to the generalized Lagrangian ([1]) with arbitrary V(9) which have the 
following form: 

V • E = Kp • B, 

9 E-VxB = K ( R B + pxE), 

V-B = 0, U 
d B + V x E = 0, 

U9 = -kE ■ B + F. (3) 



Here B and E are vectors of the magnetic and electric fields whose components are 

\e 0abc F bC} and F = -% , □ = d 2 - 



expressed via F^ as E a = F 0a , B a = -U 0abc F bc , and F = □ = d 2 



V 2 , V a = d a = £- a , a = 1,3. 

Will search for solutions of system (T5]), ([3]) with k — 1. First we notice that 
equations ([2]), ([3]) are invariant with respect to discrete transformations of space 
reflections x a — > —x a , E a — > —E a , H a — > H a , 9 — > —9 provided F is an odd function 
of 9. In other words, 9 transforms pseudoscalar. 

We will consider also the following system 

V ■ E = Kp ■ E, 

d E — V x B = n(p E - p x B), 

V • B = 0, (4) 
d B + V x E = 0, 

U9 = k(B 2 - E 2 ) + F 

which model the generalized axion electrodynamics with a scalar axionic field. 

We will present symmetries of equations ([2]), ([3]) and 01]) with arbitrary function 
F(9) and also their exact solutions, which can be found using these symmetries. We 
shall present also some results of group classification for more general systems with 
F being an arbitrary function of 9 and p^. 
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3 Exact solutions: definitions and examples 



3.1 Algorithm for finding group solutions 

Since the system ([5}, © admits rather extended symmetries, it is possible to find a 
number of its exact solutions. The algorithm for construction of group solutions of 
partial differential equations goes back to Sophus Lie. Being applied to system 
(jHJ) it includes the following steps (compare, e.g., with [15]): 

• To find a basis of the maximal Lie algebra A m corresponding to continuous local 
symmetries of the equation. 

• To find the optimal system of subalgebras SA^ of algebra A m . In the case of 
PDE with four independent variables like system (J2J), (J3]) it is reasonable to 
restrict ourselves to three-dimensional subalgebras. Their basis elements have 
the unified form Qi = ^d^ + ipfd^, i = 1, 2, 3 where Uk are dependent variables 
(in our case we can chose u a = E a , u 3+a = B a , u 7 = 9, a — 1, 2, 3). 

• Any three-dimensional subalgebra SA^ whose basis elements satisfy the condi- 
tions 

rank^f} = rank{£\ rf} (5) 

and 

rank{£f } = 3 (6) 

gives rise to change of variables which reduce system (J2]), © to a system of 
ordinary differential equation (ODE). The new variables include all invariants 
of three parameter Lie groups corresponding to the optimal subalgebras SA^. 

• Solving if possible the obtained ODE one can generate an exact (particular) 
solution of the initial PDE. 

• Applying to this solution the general symmetry group transformation it is pos- 
sible to generate a family of exact solutions depending on additional arbitrary 
(transformation) parameters. 

The first step of the algorithm presupposed finding the maximal Lie symmetry of 
the considered systems. These symmetries are presemted in the following subsection. 

3.2 Group classification of system ([2]), ([3]) 

The system of equations fl3]) includes the arbitrary element F(6) thus its symmetries 
might be different for different F. The group classification of these equation consists 
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in complete description of their continuous symmetries together with the specification 
of all functions F corresponding to different symmetries. 

Applying the classical Lie algorithm (refer, e.g., to [IS]), to find symmetries of 
system fl2]), fl3]) w.r.t. continuous groups of transformations B — > B', E — > E', 9 — > 
9', x^ — > x'^ we introduce the infinitesimal operator 

Q = ^d, + rfd B i + C j d Ej + ad 9 , (7) 
and its prolongation 

■ d ■ d 

Q(2) Q ■ i//— ■ (,,'-— • aide. + o ik d Blk (8) 

where B\ = diB j , Ej = diE j , 0j = di9, 9 ik = d$ k . Functions rf iy a ik can be 

expressed via £ l , , (i , a using the following relations: 

^ = dm) - BiD t (e), d j = Di(c j ) - Eim h ), 

<n = Di(a) - 4A(^), *ik = D k {°i) ~ duD k {i l ) 

where A = 9s + B{d & + E\d E3 + 9A + 9 ik de k . 

The invariance condition for system (jSJ), ([3]) can be written in the following form: 

Q{2)J z W=o = o (9) 

where J 7 is the manyfold defined by relations (|2J), ([3]). Then, equating coefficients for 
linearly independent functions E\ B\ 9 and their derivatives we obtain the following 
system of determining equations for the coefficients of the infinitesimal operator: 

e Ba = o, e Ea = o, c = o, &+& = o, (io) 

&E a = 0, <7 B a = 0, a ee = 0, (11) 
Da + (a - 2£° ) (F + kE a B a ) - K(B a ( a + E a r] a ) - aF e = 0, (12) 
- 2^ = 0, (13) 



— £abcV E b = 0; £,x° ~ e abcC B b = 0, 

daV a = 0, 9aC a + B a d a a = 0, 

t& + e abc Q b = 0, Co + B a a x o - e abc {v c x * + E b a x c) = 0, (14) 
V a + B a a e + Q - B\ a Eb + e abc E b e xC = 0, 

Vb« - rj b Bb = 0, rj% a - (e* = °> te ~ B a rf Eb = 0, Q - E a v b Eb = 0. 
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Here the subscripts denote the derivatives with respect to the corresponding variables: 
£g a = g|r£, etc, and there are no sums over repeating indices in the last line of ffT4|) . 

In accordance with equations ffTOj) functions £ M are Killing vectors in the space of 
independent variables: 

f = 2x»fx v - j»x v x v + c tMU x v + dx» + (15) 

where / M , d, e M and = — c Ufl are arbitrary constants. 

In accordance with ffTTj) a = ^l^ 1 + <£>2> where (p± and y?2 are functions of x^. 
Substituting this expression into f[T2l we obtain: 

y^F, + <p 2 Fg + 2 (C° -<pi)F + 2k (C° - ^) £ a 5 a 

+ K (5 a c a + - en?! - n^ 2 - 2^8^ = o. 1 j 

If the terms 8Fg, Fg, F, and 1 are linearly independent, then it follows from 
tfH| that 

a = C = 0, ITC a + £V = 0. (17) 

Substituting (ITT]) and (JTBT) into ({TBI we obtain that / w = 0, hence (TT5]) reduces to the 
form 

f = cTx v + (18) 

It follows from (fill) . (ITTjl and <IIB]) that 

^ = c^^ 6 + e abc c ob E c , ( a = c ab E b - e abc c ob B c . (19) 

In accordance with (fT8l) and (fl9i) the infinitesimal operator (JJJ) is reduced to a 
linear combination of the following basis elements: 

-Fb — 8q, Pa = 8 a , 

J ab = x a d b - x b d a + B a 8 Bb - B b d Ba + E a 8 & - E b d Ea , (20) 
Joa = x d a + x a d + e abc (E b d B c - B b 8 E .) 

where e abc is the unit antisymmetric tensor, a,b,c— 1, 2, 3. 

For some particular functions F, namely, for F — 0, F = c and F = bexp(a6) 
the symmetry of system (T5]), Q is more extended. The corresponding Lie algebra 
includes the following additional basis elements: 

p 4 = dg, D = x d + xA - B i d Bl - E { d & if F{6) = 0, (21) 
P 4 = 8 e if F{9) = c, (22) 
X = aD- 2P 4 if F{6) = be ad . (23) 
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Using the standard algorithm of group classification (see, e.g., [T5]) we can find 
symmetries of a more general system (|2J), (EJ) with arbitrary element F being a func- 
tion of both 9 and its derivatives p M . Restricting ourselves to the case of Poincare- 
invariant systems we find that F can be an arbitrary function of 9 and p^ . Moreover, 
all cases when this symmetry can be extended are presented by the following formulae: 



where /(.) is an arbitrary function on the argument given in the brackets and k is an 
arbitrary constant. Symmetry algebras of system ([2]), ([3]) where F is a function given 
by formulae ([21]), fl25J) and (1251) include all generators (j2DJ) and operators presented 
in (12Tj) . (1221) and (1251) correspondingly. 

Finally, the group classification of equations (jlj) gives the same results: this system 
is invariant w.r.t. Poincare group for arbitrary F . System (jlj) admits more extended 
symmetry in the cases enumerated in equations (|2Tj) ([26|) . 

3.3 Optimal subalgebras 

Thus, to generate exact solutions of system (j2J), (jHJ) we can exploit its invariance 
w.r.t. the Poincare group whose generators are presented in equation (TSOj) . The 
subalgebras of algebra p(l,3) defined up to the group of internal automorphism has 
been found for the first time in paper [T7]. We use a more advanced classification 
of these subalgebras proposed in |18j . In accordance with [T8] there exist 30 non- 
equivalent three-dimensional subalgebras Ai, A2, ■ ■ ■ A 3 q of algebra p(l,3) which we 
present in the following formulae by specifying their basis elements : 



F = K-Pfj.p' 1 , 

F = e ae f (p^e- ae ) 



(24) 
(25) 
(26) 



A 4 : 

A 7 : 

A 9 

A u 

A 13 

A 15 

A 17 

A 19 

A 2 i 

A22 

A 2 4 

A 2Q 

A 2 * 



Po,P u P*)\ A 2 : (Pi,P 2 ,P 3 ); ^3: (P -P 3 ,Pi,P2); 

J 03 ,P 1 ,P 2 ); A 5 : (J 03 ,P -P3,Pi); A 6 : (J 03 + aP 2 , P , P 3 ); 

J 03 + aP2,P -P 3 ,Pi); A s : (J 12 ,P ,P 3 ); 

J 12 + aP ,Pi,P 2 ); A w : (J 12 + aP 3 ,Pi,P 2 ); 

{J 12 — Po + P3, Pi, P2); A12 : (Gi, Po — P 3 , P2); 

(G 1 ,P -P 3 ,Pi + aP 2 ); A u : (G x + P 2 , P - P 3 , Pi); 

(d - P , P - P 3 , P 2 ); A 16 : (Gx + P , Pi + aP 2 , P - P 3 ); 

(J 03 + aJ 12 ,P ,P 3 ); A 18 : (aJ 03 + J 12 , Pi, P 2 ); 



(27) 




^12 — Po + P3; Gi, G 2 )\ A 2 7 : (J03 + 0iJ\2, Gi, G2); 



Gi, G2, J12)', A29 : (J01, J02, J12)] A 30 : (Ji 2 , J 2 3, J31). 



6 



Here P^ and J^ u are generators given by relations f )2"0"]) . G\ = J i — J13, G 2 = J02 — 
J 23 , a and (3 are arbitrary parameters. 

Using subalgebras ( 127|) we can deduce exact solutions for system ([2]), ([3]). Notice 
that to make an effective reduction using the Lie algorithm, we can use only such 
subalgebras whose basis elements satisfy conditions ([5]). This condition is satisfied 
by basis element of algebras A\ — A 27 but is not satisfied by A 28 , A 29 , A 30 and A e 
with a = 0. Nevertheless, the latter symmetries also can be used to generate exact 
solutions in frames of the weak transversality approach discussed in [19] . We will use 
also a certain generalization of this approach. 

In the following sections we present the complete list of reductions and find exact 
solutions for system (T2]), ([3]) which can be obtained using reduction w.r.t. the sub- 
groups of Poincare group. We will find also some solutions whose existence is caused 
by symmetry of this system with respect to the extended Poincare group. 



3.4 Plane wave solutions 

Let us find solutions of system ([2]), 03]) which are invariant w.r.t. subalgebras A\, A 2 
and A 3 . Basis elements of all these subalgebras can be represented in the following 
unified form 

A : (Pi,P 2 ,kP + sP 3 ) (28) 



where e and k are parameters. Indeed, setting in ( 128)1 e = — k we come to algebra 
A3, for e 2 < k 2 or k 2 < e 2 algebra (128)1 is equivalent to A\ or A 2 correspondingly. 

Starting from this point we mark the components of vectors B and E by subindices, 
i.e., as B = (B u B 2 , B 3 ) and E = (E u E 2 , E 3 ). 

To find the related invariant solutions we need invariants of the group whose gen- 
erators are given in (128)1 . These invariants include the dependent variables E a ,B a ,6 
(a = 1,2,3) and one independent variable u = exq — kx 3 . Let us search for solutions 
of 02]), 02]) which are functions of u. Then equations 02]) are reduced to the following 
system: 

B 3 = 0, E 3 = 6B 3 , kE 2 = -eBx, kE x = eB 2 , 

eE-y - kB 2 = 6(kE 2 + eBj), kB x + eE 2 = 6(eB 2 - kE x ) 1 ' 

where B 3 = ^f. 

The system Q29l) is easily integrated. If e 2 = k 2 7^ then 

E, = jB 2 = Fx, E 2 = -jB x = F 2 , E 3 = c9 + b, B 3 = c (30) 
k k 

where F\ and F 2 are arbitrary functions of u while c and b are arbitrary real numbers. 
The corresponding equation (J3J) is reduced to the form e 2 6 = F{6) — be, i.e., 6 is 
proportional to F(8) — be if e 7^ 0. If both e and F equal to zero then 9 is an 
arbitrary function of u. 
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For e 2 ^ k 2 solutions of f[2"9l have the following form: 



Bi = kc\Q — kbi + ec 2 , B 2 = kc 2 9 — kb 2 — eci, B 3 = c 3 , 
Ei = ec 2 9 — eb 2 — kci, E 2 = —ec\0 + eb\ — kc 2 , E 3 = c 3 9 — b 3 (e z — k 



2 ^ (31) 



where b a and c a (a = 1, 2, 3) are arbitrary constants. The corresponding equation ([3]) 
takes the form 

= -b6 + c + F (32) 

where b = (c\ + cl + ^zpj > ^ = (i^rpy and c = Ci&i + c 2 b 2 + c 3 b 3 . 
If F = or F = —m 2 8 then (132]) is reduced to the linear equation: 

= -qB + c (33) 
where a = c 2 + c 2 + flzp • Thus 

9 = a M cos /iw + 6^ sin jiuj H — -, if a = /x 2 > 0, (34) 



/x 2 



<p = a a e™ + 6 CT e _aw - — if a = -v 2 < 0, (35) 



= 2 CCj2 + ClW + C2 ^ a = (36) 

where a M , b^, a a , b a , c\ and c 2 are arbitrary constants. 

One more plane wave solution of equations (J2J), © with « = 1 and F = can be 
written as: 

E 2 = Cue cos u + sin u, E 3 = c^e sin u — d^e cos u, 

B 2 = Ckksmuj — d^k cos u, B 3 = — c^k cosu — d^k sin u, (37) 
Ex = e, B\ = 0, 6* = axo + vx\ + c 3 , uj = sxq + kx\ 

where e, d^, 5, fc, a, i/ are arbitrary constants restricted by the only relation: 

e 2 - k 2 = ve - ak. (38) 



If e = k then a = v and formulae (|37j) present solutions depending on one light 
cone variable Xq+x\. However, for e ^ k we have solutions depending on two different 
plane wave variables, i.e., ex + kx\ and ax + vx\. 

It is interesting to note that for fixed parameters a and v solutions (|37|) for E a and 
£> a satisfy the superposition principle, i.e., a sum of solutions with different e,k,Ck 
and dk is also a solution of equations (EJ, (EJ) with k = 1 and F = 0. 

Using symmetries of system ([2$, (J3]) it is possible to extend the obtained solutions. 
Indeed, applying to ( 13"Uj) . (13"Tj) the rotation transformations 

Ea E' a = RabEb } B a ^ B' a = RabBb, (39) 
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where {R a b} is an arbitrary orthogonal matrix of dimension 3x3, and then the 
Lorentz transformations 



j/sinhA i \ \ 771/ 1— cosh A 



E' a -> E' a cosh A - Eabc \ b B' c ^ + X a X b E' b — - — . 

(40) 

B' a cosh A + e abc X b E' c ^ + A a A bJ B[^f^, A = yf\\ + X 2 2 + A| 



and transforming u; — >■ n^x' 4 where are components of the constant vector given 
by the following relations: 

n = cosh A - uX a R a3 Smh - , 

(41) 

sinhA (1 — cosh A) 
n a = ^-K a 3 — A a — uX a X b R b3 — , 

we obtain more general solutions of equations (J2]), (J3j). 

In formulae (l39l) - (|4Ti) summation is imposed over the repeated index b, b = 1, 2, 3. 



3.5 Selected radial and cylindric solutions 

Let us present some exact solutions of equations (j2]), ([3]) which can be interesting 
from the physical point of view. 

First we consider solutions which include the field of point charge, i.e. 

E a = q^§, a = 1,2, 3 (42) 

where r = \/ x\ + x\ + x\ and q is a coupling constant. Notice that up to scaling the 
dependent variables x a we can restrict ourselves to q — 1. The related vector B a is 
trivial, i.e., B a = 0, while for 9 there are two solutions: 

9 = ^ «n&e = ± : (<p 1 {x + r) + <p 3 (x -r)) (43) 

where ipi and ^2 are arbitrary functions of xq + r and xq — r correspondingly, c a are 
arbitrary constants and summation is imposed over the repeating indices a = 1, 2, 3. 
These solutions correspond to trivial nonlinear terms in (j2J), ([H]). 

Radial solutions which generate nontrivial terms in the r.h.s. of equations 02]), ([2]) 
with F = —m 2 9 can be found in the following form: 

D Q x a jp q9x a . . , _ 2 . . 

S Q = -, E a = — , 9 = ci sm(mx )e - (44) 

where ci and q > are arbitrary parameters. The components of magnetic field B a 
are singular at r = while E a and # are bounded for < r < 00. 

Solutions ( 14"2"]) -( 14"4|) where obtained with using invariants of algebra A 3 q. 
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Let us present solutions which depend on two spatial variables but are rather 
similar to the three dimensional Coulomb field. We denote x = y/x[+~x%, then 
functions 

Ei = -B 2 = % E 3 = 0, B X = E 2 = % B 3 = b, 9 = arctan f — ) (45) 
x 6 x 6 \ x i J 

where b is a number, solve equations (j2J), ([3]) with k — 1 and F = 0. 

A particularity of solutions ( l4"5j) is that, in spite of their cylindric nature, the 
related electric field decreases with growing of x as the field of point charge in the 
three dimensional space. 

Functions (I45p solve the standard Maxwell equations with charges and currents 
also. However, they correspond to the charge and current densities proportional to 
1/x 3 which looks rather nonphysical. In contrary, these vectors present consistent 
solutions for equations of axion electrodynamics with zero axion mass. 

Solutions f T4B]) can be expressed via invariants of the subgroup of the extended 
Poincare group whose Lie algebra is spanned on the basis (P , P 3 , J i2 + -P4), see 
equations (120]) . ( 122]) for definitions. 

Let us write one more solution of equations ([2]), (|3]) with F = 0: 

R r> X 2 X 3 R X ( X \ 

Bi = -r— , B 2 = B 3 = — -, 6 = arctan — , (46) 

r z x r z x r z \x 3 J 

E a = ^, a = 1,2,3 (47) 

where r = \J x\ + x\ + irf, x = \Jx\ + x\. The electric field (T4T]) is directed like the 
three dimensional field of point charge but its strength is proportional to 1/r instead 
ofl/r 2 . 

Let us note that functions (T461) . (T4TI) solve equations (jlj) with k = 1,F — also. 
Two more stationary exact solutions for these equations can be written as: 

E a = % a = 1, 2, 3; 5 a = 0, 9 = ln(r) (48) 



and 



E a = ^, B a = b a , 9 = ln(r) (49) 



where b a are constants satisfying the condition b\ + b\ + 6| = 1. Functions f l4"9]) solve 
equations (T4J) with F = 0for0<r<cxD while formula (T4H]) gives solutions of equation 
® with F = p a p a . 

The complete list of exact solutions for equations (T5]), (j3J) obtained using symme- 
tries w.r.t. the 3-dimensional subalgebras of the Poincare algebra is presented in the 
following section. 
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4 Complete list of invariant solutions 



In this section we present all exact solutions for equations (J2J), (J3J) which can be ob- 
tained using symmetries w.r.t. the 3-dimensional subalgebras of the Poincare algebra. 
Basis elements of these subalgebras are given by relations (|27|) . 

We shall consider equations (EJ), (EJ) with the most popular form of function F, 
i.e., F = —m 2 6, which is the standard choice in axion electrodynamics. In addition, 
up to scaling the dependent variables, we can restrict ourselves to the case k = 1. 
Under these conventions the system ([211 , can be rewritten in the following form: 

V • E = p B, 

<9 E — V x B = p B + p x E, (50) 

V ■ B = 0, <9 B + V x E = 0, 

U9 = -E-B-m 2 # (51) 

In the following we present exact solutions just for equations (150|) . (1ST]) for both 
nonzero and zero m. 

Solutions corresponding to algebras A\ — A 3 have been discussed in previous sub- 
section. Here we apply the remaining subalgebras from the list ( 12 7p . grouping them 
into classes which correspond to similar reduced equations. 



4.1 Reductions to algebraic equations 

Let us consider subalgebras A u , A 20 , A 26 and show that using their invariants the 
system ([50]) . (15T1) can be reduced to algebraic equations. 
Algebra A n : ( J 12 - P + P 3 , A, P2) 

Invariants / of the corresponding Lie group are functions of the dependent and 
independent variables involved into system ( 150]) . ( |5T|) . which satisfy the following 
conditions 

P 1 / = P 2 / = 0, ( J 12 - P + Ps)I = 0. (52) 

The system ( )52|) is non-degenerated thus there are eight invariants which we choose 
in the following form: 

Ii = Ei sin ( — E 2 cos (, I 2 = E 2 sin £ + E\ cos £, 

J 3 = Pi sin Q — B 2 cos C, -^4 = B 2 sin C + -Bi cos (, (53) 
J 5 = P 3 , I & = B 3 , I 7 = 9, I s = u = x + x 3 

where £ = |(x 3 — x ) and I a , a = 1, 2, ...7 are arbitrary functions of w. Solving (I55|) 
for P a , P a and ^ and using (1501) we obtain 

Ei = B 2 = ci sin ( + c 2 cos P2 = —B\ = c 2 sin ( — c\ cos C, (54) 
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E 3 = c 3 9 + C4, B 3 = c 3 , (55) 

where ci, C2, C3, C4 are arbitrary real constants and is a function of w which, in 
accordance with (1511) . should satisfy the following linear algebraic relation: 

(c\ + m 2 )9 + c 3 c 4 = 0. (56) 

Thus 9 = — i 3 , Ci g if the sum in brackets is nonzero and 9 is an arbitrary constant 
provided c 3 = m = 0. 

Notice that solution (15^1) can be generalized to the following form: 

E 1 = B 2 = /(C), E 2 = -B l = (7(C), £3 = c 3 9 + c 4 , H 3 = c 3 (57) 

where /(C), <?(C) are arbitrary functions and 9 again is defined by equation ( l56]h 
However, solution (j57p cannot be obtained via symmetry reduction. 

In analogous way we obtain solutions corresponding to subalgebras A 20 and A 2Q . 
Algebra A 20 : (G u G 2 , P - P 3 ) 

c 2 -2c l xix 2 + c 2 {x\ - x\) + 2c 3 xi + 2c 3 x 2 9 + 2c 4 x 2 

Bl = E 2 = ^ \~ Pi, 

uj 2u d 

_ „ ci Ci(x 2 - x\) + 2c 2 xix 2 + 2c 3 x 2 - 2c 3 Xx9 - 2c A x 1 

B 2 = —hj X H = ~ 7, V (f2, 

u 2bj 6 

-Cix 2 + c 2 xi + c 3 -C\X X - c 2 x 2 + c 3 9 + c 4 

-D3 = 5 , E 3 = 

U3 Z uj z 

where ipi are functions of uj = xo + x 3 , 

(Ci^i + C 2 <f 2 )uJ 3 + C 3 Cj 2 2 

9 = 5 7T—. if Co + m > 0; 

9 = ip 3l cxipi + c 2 ip 2 = if C3 + m 2 = 0. 
Algebra A 26 : (J 12 - P + P 3 , d, G 2 ) 

ClX!X 2 , c 2 x 2 ex ((# - 2) W 2 + 2(x 2 - x 2 ) 

Bl = C ° SC + ^ + 4^ SinC ' 

Cl x 4 x 2 <* (((9 - 2V 2 - 2(x? - ^1) 
#2 = 3— sin C H t— 3 cos C, 

£3 = — o-sinC r H 2~ cos ^ E x = -B 2 cosC, 

U) z UJ 6 UJ Z UJ 

J? D , Cl • f T? ClX2 • /• , C2 , ClXl A 

E 2 = B X A sin C, E 3 = — — sin C H n H 5- cos C, 

UJ UJ Z UJ Z U3 Z 

9 = 0, if m^O, 9 = <p(u) if m = 
where C = 77 + f , ^ 2 = x\ — x\ — x\ — x\ and <^( w ) is an arbitrary function of 

U = X + X 3 . 
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4.2 Reductions to linear ODE 

The next class includes subalgebras A 5 , A 7 , Ai 5 , Aiq and A 25 . Using them we shall 
reduce the system (1501) . (|5T|) to the only linear ordinary differential equation (|33l) . 

Let us start with algebra A 5 whose basis elements are (J 03 ,Po — P3,Pi)- The 
corresponding invariant solutions of equations (I50p . (15 ip have the following form: 

B 1 = E 2 = [x + x 3 ) (a0 + c 2 ) , B 2 = -Ex = ci(x + x 3 ), 
B 3 = -c 3 9 + c 4 , E 3 = c 3 , c x c 2 = 0. 

Function 9 = <p(uj) depends on the only variable u = x 2 and satisfies equation 
fl33|) where a = c 2 — m 2 , c = c 3 c 4 . Its possible explicit forms are given by equations 

([MD-flMD. 

Algebra A 7 : ( J 03 + aP 2 , P - P 3 , P x ) 

D -c x + c 2 -ac 3 6> + ac 4 - a 

x + x 3 x + x 3 

B 3 = -c 3 9 + c 4 , E 3 = c 3 . 



Possible functions 9 = (p(oj) again are given by equations fl34|) - fl36|) where a 

3 



Co — m 2 , c = c 3 c 4 and oj = x 2 — a In \xq + x 3 



Algebra A 15 : (G\ - P , P - P 3) P 2 > 

Bi = E 2 = -c 2 (x + x 3 )9 - ci(x + x 3 ), B 3 = c 2 9 + ci, 

B 2 = -E x = c 3 (2u - xi) + c 2 (x + x 3 ), E 3 = c 3 (x + x 3 ) + c 2 



where uj = X\ + \{xq + x 3 ) 2 . The possible 9 are given by equations ( 135|) . (136]) where 
a 2 = m 2 + c 2 , c = cic 2 . 

Algebra A l6 : (G?i + P , Pi + «P 2 , P - P 3 ) 



Pi = (zo + ^3)(c 3 6 l - c 4 ) + ^ci(x + x 3 ) 2 + - ° 5 J 9 - a) + c 2 , 

2 1 + cr 



P 2 = ci(w - ^(x + x 3 ) 2 ) + c 3 (x + x 3 ) + - ° 5 - (afl + 1) + ac 2 , 
2 1 + cr 

P3 = -c 3 - ci(x + x 3 ) + c 4 , Pi = -P 2 - aci, 

P 2 = Pi + Ci, P 3 = -aci(x + x 3 ) + c 3 

where a; = x 2 — axi — %(xq + x 3 ) 2 , 

1 f c 2 . 1 \ 
6> = -5— — -f^ 3 + -(c 3 c 4 + cic 5 )w 2 + c 7 w + c 8 if c 2 =m 2 , 
cr + 1 \ 2 / 

c 2 a; 



L 1 9/9 

9 = ¥ + 2 if c 3 ^ m 



c 3 — m 2 
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2 2 

Here (p is the function of a; given by equations where /i 2 = —a 2 = , c = 

C3C4+C1C5 

° Algebra A 25 : ( J 03 + «?i + 0P 2 , G X ,P Q - P 3 ) 

o P c 3 + M + c 2 )C „ „ /3c 3 # + Cl + c 3 C 
-t>i = &2 = ; , -t>2 = — &\ = ; , 

X 3 + X X 3 + X 

B 3 = c 3 9 + c 2 , E 3 = -c 3 , 

where ( = x\ — a In \x 3 + x | and 6 1 = (p(u>) is a function of u = x 2 — /31n |x 3 + x | 
given by equations (13"4")) - ([3T))) with c = — c 2 c 3 and /i 2 = —a 2 = c 2 — m 2 . 

Consider now reductions which can be made with using invariants of subalgebras 
At, A 8 , Aig, A 2 4 and A 27 . In this way we will reduce the system fl50l) . fl5T|) to linear 
ODEs which, however, differ from (133]) . 

Algebra A 4 : (J 03 ,Pi,P 2 >. 

-c 2 x 3 6» + c 6 x 3 - cix -cix 3 6» + c 5 x 3 + c 2 x 

-Dl = 5 , -D2 = 5 , -t>3 = C 3 , 

fl fl (58) 

-CiXofc* + C 5 X + C 2 X 3 C 2 X 9 + CiX 3 - C 6 X 

Ei = = , E 2 = , E 3 = c 3 6 + c 4 

UT UJ Z 

where c\, ■■ ■ ,Cq are arbitrary constants, 9 = 6{uj) and uj 2 = x\ — x\. Substituting 
([58]) into fl5T|) we obtain: 

u 2 6 + QU0 + {y 2 + jj 2 uj 2 )6 = 5 + au 2 (59) 



where v 2 = c 2 + c 2 , fi 2 = c 2 + m 2 , 5 = C1C5 + c 2 C6, a = c 3 C4 and # = 86 /dui. 

c 2 > x 2 



The general real solution of equation fl59|) for X 2 , > x§ is: 



= c 7 (J iu (fiuj) + J_ iv (yua;)) + c 8 (^(/iw) + Y_i„(/4u;)) 

( coth (~2~) Jiu ^ u ^ + iE ^(f iU )) + -^-M 1 ? iz/ > A^) ( 60 ) 



where a; = \/xq — x 3 , Ji u (/j,uj) and Y- w (/iuj) are Bessel functions of the first and second 
kind, L s (l, iv, fiu) is the Lommel function s, J iu (fiu) and E iu (fiu) are Anger and 
Weber functions. 

If \w = and x\ > x\ then solutions of (159]) are reduced to the following form: 



5 otuj 2 

6 = c 7 sin(z/ In w) + c 8 cos(z/ In w) H — - H — if = 0, z/ 7^ 0; (61) 

i/ 2 + 4 



6»= I aw 2 + ^ln 2 (w) + c 7 ln(w) + c 8 if fi = v = 0; (62) 

61 = c 7 J (/iw) + c 8 Y (fiu) H — - if v = 5 = 0, fi ^ 0. (63) 

A* 
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We shall not present the cumbersome general solution of equation (IqT?|) for Xq — x\ < 

2 

but restrict ourselves to the particular case when a = ^5. Then 
9 = c 7 {l iv {iMZj) + I- iv ([tij)) + c 8 (KivQiu) + K_ iv (/xo))) + — 

7 1^ 



where 00 = \Jx\ — x^. 

Algebra A 8 : {J 12 ,P Q: Pz) 

C 2 X 2 9 + CiXx - c 6 x 2 -c 2 x x 9 + CiX 2 + CqXi 
E>\ = 5 , -t>2 = n , -t>3 = + C4, 

c x xi9 + C5X1 - c 2 x 2 cix 2 6 + c 5 x 2 + c 2 x 1 

&1 = 7, , &2 = 5 , &3 = C3 

where u 2 = x\ + x\ and 9 is a solution of equation (159]) with 

z/ 2 = c 2 — c 2 ,, /U 2 = c 2 — m 2 , 5 = C1C5 + c 2 c 6 , a = c 3 c 4 . (64) 

If c 2 > c 2 . and c§ > m 2 then 6 1 is defined by relations fl60l) - fl63|) where /1, v and 5 are 
constants given in If c 2 — c 2 , = —A 2 < 0, m 2 < c 2 and a(a\ 2 + 5/i 2 ) = then 

9 = c 7 J A (/iw) + c 8 Y A (/iw) - ^ ^cot (^Y^j J\{vw) + E A (yUu;)^ + 

where J A (/iCoO and Y A (/icu) are the Bessel functions of the first and second kind, J\(iiuj) 
and E\(fiu) are Anger and Weber functions correspondingly. In addition, 

9 = c 7 u x + c s u- x - — - |L, A 2 = c 2 . - c 2 if c 2 > c 2 , c 2 = m 2 ; (65) 
61 = c 7 I a (ku;) + c 8 K x (koj) + / if m 2 - c 2 = k 2 > 0, cj > c 2 , (66) 

where 



A fc-^ 4-rv ry 

/ = -- if 5 = a-, A^O, /= - if A = 2, 5 = 0, 

f = ~ if 5 = A = 0, 
2k 



Ia(^c^) and K a (/tXo>) are the modified Bessel functions of the first and second kind. 

Solutions (16"61) are valid also for parameters 5 and A which do not satisfy condi- 
tions presented in (|67|) . The corresponding function / in (166]) can be expressed via 
the Bessel and hypergeometric functions, but we will not present these cumbersome 
expressions here. 

Algebra A l9 : (J 12 , J 3, ^0 - ^3) 

D „ e 1 {x 1 +x 2 9) c 1 (x 2 -x 1 8) 

£>\ — &1 = 7 ; w 9 , oT; ^2 = 



(X 3 + Xq)(x\ + X 2 2 y ' ^ (X 3 + Xq){x\ + x 2 2 y 
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B 3 = -c 3 9 + c 2 , E s = c 3 , 

where 9 is a function of u — \J x\ + x\ which solves equation flol?]) with v — 5 — 

0, /i 2 = c 2 — m 2 , a = c 2 c 3 . Its explicit form is given by equations (I6T|) and f[6"6"j) 
were 5 = 0. 

Algebra A 24 ■ (G u J 03 , P2) 

B l = -x 3 ip, B 2 = ; ■ , B 3 = XxCp, 

ur x + x 3 

C 2 X 3 Ci C 2 Xi 

Ex = ^ H ■ , E 2 = x ip, E 3 = — - 

OJ 6 Xo + x 3 oo 6 

where 00 = \/ x$ — x\ — icf, y> = ^(u;). Functions and 6 1 should satisfy the following 
equations: 

u<p + 3cp+ (- + ^)e = o, 9 + — - fci + -)<p + m 2 9 = 0. 

If C1C2 = then this system can be integrated in elementary or special functions: 

c 2 9 + c 3 

v = § — ; 



9 = c 4 sinh ^— J + c 5 cosh ^— J if m = 0, c 2 7^ 0, 



Cl = 


= 


= 


c 4 




1 


= 






u; 


= 


D 









c 4 + y ^c 5 + c 2 c 3 y ^j^j dx^j if m^0, c 2 7^ 

where D = D(0, m_,n, m+,f(u))) is the Heun double confluent function with 

m± = m 2 + c 2 , ± i 71 = 2(m 2 - c 2 ), /(w) = - + - . 

4 cj z — 1 

Let c 2 = 0, ci 7^ 0, then 

if = — [9 + — + m 2 9 , 

ci y oj j 

9 = csdioj) + c 4 (G 2 + G*(u)) + ic 5 (G 2 - G*(u)) 

where 

3 + ici 3 — ici 3 m 2 a; 2 



G x {u) = F 



2 ' 2 '2' 



ici m 2 u 2 4 _ lL _ in 



G 2 (lu) = F ( 1 + ici, 1 - ici; 1 + -y ; — ) u 
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F(a, b; c; x) are hypergeometric functions and the asterisk denotes the complex con- 
jugation. 

Algebra A 27 : (Jq 3 + aJi 2 , C?i, G2) 

y?i , (x + x 3 ) 2 - xj + x 2 , xix 2 
Sl = + 2(x + x 3 )u* ^ " (xo + xs)^ 

B _ ¥2 (x + x 3 ) 2 + x\ - x\ _ x x x 2 
2 X0 + X3 2(x + x 3 )u 4 ^ 4 (x + x^u 4 ^ 3 ' 

p R , ^3(^0 + ^3) „ R (fsjxo + x 3 ) 
u; 4 a; 4 

X 2 <^3 - ^1^4 D ^1^3 + ^2<^4 
E ' = ^ ' Bi = ^ ' 

9 = — (qJ^mw) + c 2 Y 1 (ma;)) if m/0, w 2 = Xq — x 2 — x 2 . — x\ > 0, 
9 = — (cili(md)) + c 2 Ki(mu)) if m ^ 0, a) 2 = — u 2 > 0, 
6* = c\ H — 77 if m = 0, 

V?i = c 2 cos (a In (x + x 3 )) + c 3 sin (a In (x + x 3 )), 
(f 2 = c 2 sin (a In (x + x 3 )) - c 3 cos (a In (x + x 3 )), 
uj 2 \ ( , ui 2 



¥>3 = c 4 sin a In + c 5 cos a In 

X + X 3 J V ^0 + ^3 

w 2 \ . / , u 2 



V?4 = c 4 cos I a In J — c 5 sin I a In J , (c 3 + c 2 ) (c 5 + c 4 ) = 0. 

V ^0 + ^3/ V x o + X 3 J 

4.3 Reductions to nonlinear ODE 

Using subalgebras A 6 , A 9 , A 10 , A i3 , A u , An and A 1S we can reduce (1501) . (IqT|) to 
systems of ordinary differential equations which however are nonlinear. 
Algebra A 6 : ( J 03 + aP 2 , P , P 3 ) , a ^ 

n 1 x 2 X 2 X 2 X 2 

B 1 = cosh v?2 sinh — , B 2 = a(p 2 cosh a<fii smh — , 

a a a a 

B 3 = -c x 9 + c 2 , 

-Ei = aipi cosh — — a0 2 sinh — , E 2 = (pi sinh — — y?2 cosh — , E 3 — C\ 
a a a a 

where 9, <pi and (p 2 are functions of 00 — x\ which satisfy the following system of 
nonlinear equations: 

a9ip 2 = a 2 (p 2 + ip 2 , ^102 - 01^2 = c 3 , (68) 

and 

9 = (m 2 - c\)9 + a(0iipi - (fi 2 (p 2 ) + c\c 2 . (69) 



17 



We could find only particular solutions of this complicated system, which corre- 
spond to some special values of arbitrary constants. First let us present solutions 
linear in u: 



u c x c 2 cjc 5 , 

9 = 5- ± , <^2 = ± > + c 5 , </?i = c 4 v?2, c 4 < 1 (70) 

a v 1 v aJl-c 2 A 



•4 

2 ^2 



if i/ 2 = m 2 — cj > 0, 



W CiC 2 - lc 5 , ^ 2 ^ n / 71 x 

= ~ + — 2~±- , ¥?2 = ± 7= + C 5 , C 4 > 1 71 

a /i 2 fi a^c\ - 1 

if /j, 2 = cf — m 2 > 0. 

If c 2 = m 2 and v?i = ^2 then 6* is given by equation ( 136}) with c = — C1C2 while (^2 
is a linear combination of Airy functions: 

<p 2 = c 7 Ai (A(w - 1/)) + c 8 Bi (A(w - 1/)) (72) 
where A = (^) 3 , ^ = - J -. 

V « / ' aciC2 

If c 2 = m 2 , C2 = 0, C3 7^ then we find a particular solution: 

9 = ( afi 2 H — J xi + C5, <y?2 = C6 cosh /iXi + C7 sinh /iXi, <£>i = — 2 



where /1 = -^rz^s and c 2 7^ c\. 

If c 2 = m 2 and <p± = c^ip 2 , C4 7^ 1, c 2 = then 



9 = aX (pldoo + —00 + c 6 (73) 

J OL 

where A = |(c 2 — 1) and (p 2 is an elliptic function which solves the equation 

(p 2 = X(fl - mp 2 (74) 

where k = ^ z § s -. In addition to elliptic functions, equation (j74p admits particular 
solutions in elementary ones: 

V9 2 = ±W ^tanh ( W^w + c 7 j if c 5 < 1, c\ > 1, (75) 



ip 2 = ±y^tan f \J-^~u + c 7 ] if c 5 > 1, C4 < 1, (76) 

V?2 = ±-t5- if c 5 = 1, c 4 > 1. (77) 
V Aw 
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If cf = m 2 + 4y an d fi — ±a/1 + c 2 <£> 2 , then we can set 
9 = ac4f2 + CiC 2 o. 2 
and </?2 should satisfy the following equation 

1 

V?2 - C 4 02^2 + -^V?2 = 0. 

cr 

Its solutions can be found in the implicit form: 

2 r 2 dt 

LO = C 4 a / 7 ■ r h Cq 



iy ( c 2 e H^ 2 



+ 1 



where W is the Lambert function, i.e., the analytical at y — solution of equation 
Finally, for cf 7^ m 2 and <pi = <£> 2 we find the following solutions: 

■iv = yi 



9 — - (2h>ca sinh 2uu + C6 cosh 2vu;) — > 2v = ym 2 — cf , 



<£> 2 = D (0, m + , n, m_, tanh z/cu) I c 7 + c 8 



D 2 (0, m + , n, m_, tanhi/a;) 



where D(0, m + , n, m_, tanh z/w) is the Heun double confluent functions with m± = 
a ± n — If in (175]) c 5 = |^ and ^ = — ^ < then the corresponding 

expression for if 2 is reduced to the following form: 



cp 2 = c 7 J±_ Ue VUJ ) + c 8 Y±_ Ue™) . (79) 

VOL VOL 

Algebra A 9 : ( J 23 + aP , ^2,^3), a ^ 

£/i = Ci^ + c 2 , £/ 2 = <£i cos (f 2 sin — , £/ 3 = (^1 sin h y? 2 cos — , 

a a a a 

n n . x . . . x n ■ x ■ x 

Bi = Ci, B 2 — —capi cos h a<£> 2 sin — , B 3 = —a<pi sin a<£> 2 cos — 

a a a a 

where <pi, <p 2 and 9 are functions of u = X\ which satisfy equations ( 168]) and the 
following equation: 

9 = (cf + m 2 )9 - a(ipx(pt + ^2^2) + cic 2 . 

A particular solution of this system is <pi = c 4 y? 2 and 9, (p 2 given by equation (ITT]) 
where /1 2 = m 2 + cf. 

If c 2 + m 2 = then we obtain solutions given by equations (175]) . (|71j) . ([76]) (1771) 
where A = — |(c| + 1), and the following solutions: 

ipi = Cq cos/iw + C7 sin/xo;, <^2 = C6 sin/iw — C7 cos/iw, , , 

= (i - a/, 2 ) W + c 5 (80) 
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where /x = 

Algebra A w : (J 12 + aP 3 , P 1 ,P 2 ),a^O 

x 3 x 3 X 3 X 3 

£>! = (^1 COS if 2 Sin , ±5 2 = V^l Sm 1~ ¥>2 cos , -D3 = Ci, 

3/ 2 3 3 3 /i 

£/i = a0i cos aif>2 sm — , E 2 = a<fii sm h a(p 2 cos — , E 3 = c\9 + c 2 

a a a a 

where (pi, <p 2 and 9 are functions of u = x satisfying (1681) and the following equation: 

9 = —(cl + m 2 )9 - a((fi<fi + 2 y? 2 ) - c\c 2 . 

If c\ + m 2 = we again obtain solutions (1501) and solutions given by equations (175|) . 
dZ3), (USD dZZD where A = -|(cl + 1). 

Algebra A 17 : (J 03 + aJ 12 , Pq, P 3 ), a/0 

B 1 = (aLp' 2 x 2 - Lp 2 xi)e~*~ w + {x x <Px + aLp' 1 x 2 )e^ UJ , 

B 2 = —(aip' 2 xi + p 2 x 2 )e~*~ w - (ct(p[x 1 - y?i£ 2 )e«~ w , B 3 = -c x 9 + c 2 , 

E x = (aip' 2 xi + Lp 2 X2)e~«~ U] - (xiaip[ - tpix 2 )e*~ u ', 

E 2 = {pLip' 2 x 2 - tp 2 X])e~ <*~ w - (a(p[x 2 + Lpxx^e^ \ E 3 = a, 

where u = | \n(xf + x 2 ), ( = arctan^- . Functions <p x = (pi(w), y? 2 = (p 2 (uj) and 
9 = 9{u) should satisfy (168)) and the following equation: 

e- 2uj 9 + (m 2 - cj)9 + 2a{^ 2 + cp^) + Cl c 2 = 0. (81) 

This rather complicated system has the following particular solutions for c\ = ±m: 

9 = -u + c 4 , (ft = c 5 , tp 2 = c 6 ; (82) 
a 



9 = I - + ak 2 ) u + c 4 , (pi = c 5 e KaJ , <p 2 = c 6 e KU) if 2c 5 c 6 k + c 3 = 0; 



9 = --cic 2 e 2ul + c A u + c 5 , tpi = 0, (p 2 = c 6 J M (fee") + c 7 Y M (A;e u 



1 / C\C 2 2 

u = —y/c^a — 1 if = k > 0, 

a 2a 

and 

ipi = Kip 2 , 9 = </' 2 e aJ ip 2 = 2/xtan (/ze^) + C4 (83) 

if n = cic 2 = 4/z 2 > 0, a = ±1. In ( 183|) we restrict ourselves to the particular 
value of a in order to obtain the most compact expressions for exact solutions. 
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An exact solution of equation (|68[) . (j5T|) for m 2 — cf= AX 2 > and C2 = is given 
by the following equation: 

e 4^t(l+a 2 )o;+2A 2 e 2 " 

9= JeM^ + 2X^ di0 + c ^ <P2 = 0< Vi=W* (84) 
where A, fi and a are arbitrary real numbers. 
Algebra A 18 : (aJ 03 + J 12) Pi, P 2 ), a^O 

B 1 = e~ 2w (((ptXo - aip 2 x 3 ) cos ( - (y? 2 £o + atpix 3 ) sin C) , 

B 2 = e~ 2uJ {{x ipi - aLp 2 x 3 ) sin ( + (x ip 2 + aip 1 x 3 ) cos C) , B 3 = c u 

Ei = e~ 2tJ ((-a<p 2 xo + y?ix 3 ) sin ( + (ay?ix + (p 2 x 3 ) cos C) , 

E 2 = e~ 2w ((aLp 2 Xo - </?ix 3 ) cos C + (acptXo + ^2^3) sin C) , E 3 = c x Q - c 2 

where lo — | ln(xQ — a£ = In (a; + x 3 ) — ln(a;o — £3), ^2 and 6 1 are functions 

of uo which should solve the system including ( |68|) and the following equation: 

e~ 2uj 6 = -(m 2 + c\)9 + a(y?i<^i + ^2^2) + cic 2 . 

Particular solutions of this system for m 2 + c 2 = are: 

0=( a&; 2 ] W + c 4 , <pi = c 5 sin(fcu;), y? 2 = c 6 cos(£xj), /t 



a / c 5 c 6 

In addition, the solutions fl82|) . f )83|) and are valid where u; — >■ ln( 

4.4 Reductions to PDE 

Finally, let us make reductions of system (1501) . (15T1) using the remaining subalgebras, 
i.e., ^6 with a = and A 2 g — A 3 q. Basis elements of these algebras do not satisfy con- 
dition ([5]) and so it is not possible to use the classical symmetry reduction approach. 
However, to make the reductions we can impose additional conditions on dependent 
variables which force equations §5§ to be satisfied. 

This idea is used in the weak transversality approach discussed in [19] . Moreover, 
in this approach the condition ([5]) by itself is used to find algebraic conditions for 
elements of matrices (p*. 

We use even more week conditions which we call extra weak transversality. In 
other words we also look for additional constraints to solutions of equations f loTij) . f[5"Tj) 
which force condition ([5]) to be satisfied. But instead of the direct use of algebraic 
condition ([5]) we also take into account their differential consequences. As a result we 
make all reductions for the system (1501) . (|5T|) which can be obtained in frames of the 
weak transversality approach and also some additional reductions. 
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Let us start with algebra A 6 with a = 0. The set of the related basis elements 
(P , P 3 , J 03 ) does not satisfies condition ([5]). If we consider this condition as an 
additional algebraic constraint to solutions of equations fl50|) . f l5T]) then components 
of vectors E and B should satisfy the following relations: 

E 1 = E 2 = B 1 = B 2 = (85) 



Substituting (joojl into fl50|) . fjoTj) and supposing that E 3 , B 3 and 6 1 depend only 
on invariants x\ and x 2 we can find the corresponding exact solutions. However, we 
will obtain more general solutions using the following observation. 

To force condition fl5]) to be satisfied it is possible to apply additional conditions 
which are weaker than (1851) . In particular, condition ([5]) should be satisfied if we 
impose a constraint on dependent variables which nullifies the term Ex0b 2 — E 2 8b 1 — 
Bi8e 2 + B 2 8e 1 in J03. Up to Lorentz transformations such constraints are exhausted 
by the following ones: 

Ex = 0, E 2 = 0; (86) 
Ex = 0, B 2 = 0; (87) 
Ex — Bx, E 2 = B 2 . (88) 

Let relations (|86|) are fulfilled and B 2 , B 3 , Ex, E 3 , 6 depend on the invariant 
variables xx and x 2 . Then the system (1501 . (|5T|) is solved by the following vectors: 

E 1 = E 2 = 0, E 3 = c 3 , Bx = B 2 = B 3 = -c 3 6 + c 2 (89) 

OX 2 OXx 

provided function 6 = 9{xx,x 2 ) satisfies the following equation 

d 2 6 d 2 6 9 0x „ , . 

^2+^2 =(™ < 2 ~ 4)6 + Cl c 2 . (90) 

and = 0(xi,x 2 ) solves the two-dimension Laplace equation: 

dx\ dx\ 
A particular solution of equation ( 190]) is: 

= X{xx)Y{x 2 ) + ^ 2 if c\^m 2 , 

c-x m ^g2) 

9 = X(xx)Y(x 2 ) + ^(xl + x 2 2 ) it c\ = m 2 

where 

X(xi) = cs^e*"* 1 + c^e"^ 1 , F(x 2 ) = c 5i(U cosfa^) + o^sin^a^). (93) 
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Here A; 2 = m 2 + fi 2 , = c\ + fi 2 , and c Sj/J , s = 3, 4, 5, 6 are arbitrary constants. The 
general solution of equation f[9"Uj) can be expressed as a sum (integral) of functions 
( 192]) over all possible values of \i and c SjA1 . 

Solutions (189]) include an arbitrary harmonic function 0. Only a very particular 
case of this solution corresponding to = Const can be obtained in frames of the 
standard weak transversality approach discussed in |19j . 

Analogously, imposing condition ( 18 7p we obtain the following solutions: 

E x = 0, E 3 = c u B 2 = 0, B 3 = - Cl 9 + c 2 (94) 

and 

E 2 = c 3 e» K ' + c§e^ ° 4 > } 

Bl = C3e ^(c4e^-c s e-^2) _ c6e i( C5 e-^ 2 _ C4 e^2) ) 

9 = s 1 (c 4 e^ 2 + c 5 e-^' 2 ) + c 7 e^ 2 + Cge"^' 2 + 



c? — m 2 



if c\ - m 2 = n 2 > 0; 

_ c ^ g^( c 4 cos VX2 — C5 sin _|_ c ^g-i(c4 cos^Z2-c 5 sin/xX2) 

_ C3e ^(c4COS!/X2-c 5 sin/xa;2) _ ^g-^ (c 4 cos i^x 2 -c 5 sin/xx 2 ) 

C ' i \ | . QjC 7 (95) 
V = X\ ( C 4 Sin UX2 + C5 COS U1X2) + Sin UX2 + C9 COS U1X2 H — 5 ;r 

^ — m" 1 



if C g - m 2 = -z/ 2 < 0; 
£7 2 = c 3 sinh ^c 4 X2 + c 5 x 2 ^ + c 6 cosh ^c 4 x 2 + c 5 x 2 

£?i = c 6 sinh ^c 4 X2 + c 5 x 2 ^j + c 3 cosh ^c 4 x 2 + c 5 x 2 

9 = xi(c 4 x 2 + c 5 ) - ^cic 2 x 2 + C7X2 + c 8 if c 2 = m 2 . 



If conditions (1881) are imposed then one obtains the solutions 

E a = B a = a = 1, 2, E 3 = c u B 3 = -c x Q + c 2 (96) 
where is a function satisfying (19T1) . and 9 is a solution of the following equation: 

+ = (m 2 - c\)9 + Cl c 2 + (c\0) 2 + («9 2 0) 2 . (97) 

Another solution corresponding to (188]) is given by equations (185]) for B\, B 2 , Ex, E% 
and E 3 = ci, _B 3 = — C]0 + c 2 , while 6* is given by equation (|9"2"]) . 
Algebra A 28 : (Gi, G 2 , ^12) 

B 1 = E 2 = t p — tt7 (cix x - x 2 (ci9 + c 2 )) , 

(x + x 3 y 
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B 2 = -Ei = - — -| — — (dx 2 + X\{c x Q + c 2 )) , 

\Xq "t 3^3 J 

a ?> — t~ — , „ \ o ? -^3 — — — , y 



(x + x 3 ) 2 ' " (x + x 3 ) 2 ' X + Xi 



9 9 9 9 



where </2 is a function of two variables u = ' 2 (x +x 3 ) 3 an< ^ £ = ^o + ^s, which satisfies 
the following equation: 

<9V c\ 2 \ CiC 2 

-i -m )ifi+ — . (98) 



ft^C VC 4 r C 

Let m = Ci = then = y?i(w) + y?2(C) where <y?i and <£> 2 are arbitrary functions. For 
Ci = 0,m 2 7^ equation (1981) admits solutions in separated variables: 

if = ^(a M sin(z/ M £ + ) sin(/i£_) + 6 M cos(z/ M £ + ) cos(/x£_) 

+c M cos(z/ Ai ^+) sin(/j£_) + d„ sin(z/ M £ + ) cos(/x£_)) (99) 

where £± = w ± £, = m 2 + /i 2 and /x, S^, a M , 6 M , c M and are arbitrary constants. 
For ci ^ we obtain: 

Algebra A 29 : ( J 01 , J 2, ^12) 



x 2 (ci6 + c 2 ) Xi(ci9 + c 2 ) ci^o 



-Bi = 5 > ^2 = 5 > B 



C X X 2 C\X\ X (ci6 l + C2) 

-fri = 5-, -£^2 = — 5-, -£^3 = 5 , y = — 

w 3 ur or a; 

where u 2 = x^ — x 2 — X 2 . and is a function of lo and x 3 which satisfy the following 
equation: 

a5-a^ = V^ + m J v, + ^ lf *°>^ + * 2 (100) 



where u = y Xq — x 2 — x 2 , and 

S + S=(™ 2 -|)"-^ -4<4 + 4 (101) 

where a) = a/x 2 + i| — a; 2 ,. 

Let ci = m = then the general solution of equation fllOOH is: (p = <fi(co + x 3 ) + 
(p 2 (oj — X3) where ipi and ip 2 are arbitrary functions. Solutions which correspond to 
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Ci = 0, m 7^ can be obtained from by changing £ + — y x 3 , — > u. If ci ^ 
and m^fl then 

V 2 = J^ ^ (( a ^ + b ^ j sin (/^3) + + rf^y ^-J cos(/ix 2 



- ClC2 /(^/^l^ (102) 



where D M = D(0, /c^ , s, /(a;)) is the double confluent Heun function with fe- 
rn 2 — fx 2 + cf ± |, s = 2(m 2 — /U 2 — c 2 ), /(a;) = j^zj, A 4 , fyu, c M and d M are arbitrary 
constants. 

Solutions of equation f llOip also can be represented in the form ( 11021) where u — y u 
and 

kt = ^-m 2 + cl±- s = 2^-m 2 -cll f(u>) = %±±. 
^ 4 ur — 1 

Algebra A 30 : ( J 12 , J 23 , J 31 ) 

t-j cix a (ci6 , -c 2 )x 

-Da = = o , P = — , 



where r = y x 2 + x\ + x| and <^ is a function of r and xq satisfying the following 
equation: 

d 2 <p d 2 (p f cj 2 \ Cl c 2 

Solutions of this equation can be represented in the form (11021) where u = r, 
X3 — y Xo and D /t = D(0, k~, s, k+ f(u))) is the double confluent Heun function with 

k ± = _ (m 2 + ^ + c?) ± 1 j s = 2(c 2 _ m 2 _ ^ /(w) = /(r) = 

A special solution of equation (11031) corresponding to c 2 = and zero constant of 
variable separation is given in 



4.5 Solutions with maximal arbitrariness 

Solutions considered in the above include arbitrary parameters and in some cases 
even arbitrary functions. At the end of our analysis a special class of solutions will be 
presented which depend on six (!) arbitrary functions. This class cover all reductions 
which can be obtained using subalgebras Ai 2 — A u and A 2Q — A 23 . 
Let us define 

B x = E 2 = if>i(x u x 2 , u) - xi(fi(u) - x 2 (<p4,(u) - (p 5 (u)0 \ , 

B 2 = -Ei = i) 2 (x 1 ,x 2 , u) - x 2 (p 5 (u) + xi (ip 2 {u) ~ > ^ 104 ' ) 
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where ipi,... <^5 and ipiifa are functions of on = x + x 3 and xx,x 2 ,u> respectively, and 

= , if F — —m 0, m f: 0, 

m 2 

= <p 3 {u), (<Pi + <Pb)(<P2 + <Pa) = if F = 0. (105) 

Up to restriction present in fl 105[) functions tpi,<f2 and <^ 3 are arbitrary while ipi 
and ^2 should satisfy the Caushy-Rieman condition with respect to variables Xi and 
x 2 : 

8x^ + 82^2 = 0, 8x^2-82^ = 0. (106) 

It is easily to verify that Ansatz (j!04p . f)105p does solves equations (IBTil) . (|5T|) . 
Thus any solution of the Caushy-Rieman equations f)106p depending on parameter u 
and three arbitrary functions (fx(u)), <Pi{(jS), ^3^) satisfying (11051) give rise to solution 
J323D of system (|50|), (|5B. 



5 Discussion 

The main aim of the present paper is to find families of exact solutions of field 
equations of axion electrodynamics using invariants of three parameter subgroups 
of the Poincare group. The complete list of reductions which can be made using 
these invariants together with obtained solutions are presented in sections 3 and 4. 
Among them there are solutions including sets of arbitrary parameters and arbitrary 
functions as well. In addition, it is possible to generate more extended families of 
exact solutions applying the inhomogeneous Lorentz transformations to the found 
ones. 

For such subalgebras whose basis elements do not satisfy the transversality con- 
dition (p2) we apply the week and "extra weak" transversality approach, see section 
4.4. As a result we find solutions fl89]) - fl94"|) which cannot be found applying standard 
weak transversality conditions discussed in [1~9] . 

Making reductions of equations (|2J), (j3J) we restrict ourselves to functions F linear 
in 9. However, these reductions do not depend of the choice of F\ to obtain reduced 
equations with F arbitrary it is sufficient simple to change m 2 8 — > —F{9) or even 
m 2 6 — » —F{9,p lJl p^) everywhere. 

Except a particular example given by relations (l4"6"]) - (j4l?]) we did not present exact 
solutions for equations (j3J. Let us note that reductions of these equations can be made 
in a very straightforward way. Indeed, making the gauge transformation E a — > e e E a 
and B a — > e d B a we can reduce these equations to a system including the Maxwell 
equation for the electromagnetic field in vacua and the following equation: 

U9 = k(B 2 - E 2 )e- 2e + F. (107) 
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Since reductions of the free Maxwell equations with using three-dimension subalgebras 
of p(l,3) have been done in paper [20J, to find the related exact solutions for system 
(jl]) it is sufficient to solve equation ( I107|) with B and E being exact solutions found 
in [20]. 

Solutions presented in sections 3 and 4 can have various useful applications. In- 
deed, the significance of exact solutions, even particular ones, can be rather high. 
First they present a certain information about particular properties of the model. 
Secondly, they can solve an important particular boundary value problem, a famous 
example of this kind is the Barenblat solution for the diffusion equation [21J. In 
addition, the particular exact solutions can be used to test the accuracy of various 
approximate approaches. 

The found solutions, especially those which include arbitrary functions or, like 
( 137)) . satisfy the superposition principle, are good candidates to applications in var- 
ious initial and boundary value problems of axion electrodynamics. Some of these 
solutions, e.g., (j3I) . f )34|) with \i = 1 and c\ + c| > 1, describe the wave propagation 
with the group velocity higher then the velocity of light. Moreover, these solutions 
are smooth and bounded functions which correspond to positive definite and bounded 
energy density [14] . 

We believe that the list of exact solutions presented in sections 3 and 4 can find 
other interesting applications. In particular, solutions, which correspond to algebras 
A 9 , Ax, A 17 , Ais and A 2 s generate well visible dynamical contributions to the axion 
mass. In addition, as it was indicated in [5], the vectors of the electric and magnetic 
fields described by relations fj4"5]) give rise to exactly solvable Dirac equation for a 
charged particle anomalously interacting with these fields. We plane to present a 
detailed analysis of the obtained solutions elsewhere. 
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